Extraclassical receptive field phenomena & short-range connectivity in VI 
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Abstract 

Neural mechanisms of extraclassical receptive field 
phenomena in VI are commonly assumed to result 
from long-range lateral connections and/or extras- 
triate feedback. We address two such phenomena: 
surround suppression and contrast dependent recep- 
tive field size. We present rigorous computational 
support for the hypothesis that the phenomena 
largely result from local short-range (< 0.5 mm) 
cortical connections and LGN input. Surround 
suppression in our simulations results from (A) 
direct cortical inhibition or (B) suppression of re- 
current cortical excitation, or (C) action of both 
these mechanisms simultaneously. Mechanisms B 
and C are substantially more prevalent than A. We 
observe an average growth in the range of spatial 
summation of excitatory and inhibitory synaptic 
inputs for low contrast. However, we find this is 
neither sufficient nor necessary to explain contrast 
dependent receptive field size, which usually involves 
additional changes in the relative gain of these inputs. 

keywords: visual cortex, spatial summation, surround 
suppression, receptive field, model, simulation 



Introduction 

In mammals, the very first stage of cortical visual processing 
takes place in the striate cortex (area VI). Already at this 
level spatial summation of visual input displays a considerable 
complexity. This is manifest from the fact that single cells in 
VI display surround suppression (suppression for increasing 
stimulus size, "size tuning") of their responses and have re- 
ceptive field sizes (preferred sizes) that depend on contrast 
and context. Such behavior is seen throughout striate cortex, 
including all cell types in all layers and at all eccentricities 
1-9 . The suppression seen is substantial, 30-40% on average 
in macaque VI s . Similarly profound is the growth of receptive 



field size at low contrasts. Typical is a doubling in receptive 
field size for stimulus contrasts decreasing by a factor of 2-3 
on the linear part of the contrast response function 4 . Ap- 
parently, neurons in VI sacrifice spatial sensitivity in return 
for a gain in contrast sensitivity at low contrasts 4 . Neu- 
ral mechanisms responsible for these two so called extraclas- 
sical receptive field phenomena are very poorly understood. 
Understanding these mechanisms is potentially important for 
developing a theoretical model of early signal integration and 
neural encoding of visual features in VI. 

Popular working hypotheses are that extraclassical recep- 
tive field phenomena in VI are a product of long-range hori- 
zontal connections within VI 10-13 and/or feedback from ex- 
trastriate areas 6 . 8 . 14 > 15 . Arguments in support of these hy- 
potheses are based on the observed surround sizes and the cor- 
tical magnification factor, and claim that short-range (< 0.5 
mm) and even long-range horizontal (< 5 mm) connections 
in VI do not have sufficient spatial extent to be responsible 
for surround suppression or receptive field growth 6,s . Fur- 
ther support along this line was presented using anterograde 
and retrograde tracer injections 1 and timing experiments 15 . 
So far, however, all support for these hypotheses is entirely 
based on indirect experimental observations, while they also 
lack any rigorous computational support. 

The hypothesis that the phenomena result from local short- 
range (< 0.5 mm) cortical connections and LGN input is 
largely ignored or dismissed. However, support for it can 
be found in the experimental data. For instance, surround 
suppression and contrast dependent receptive field size are 
equally profound throughout VI 4,6 , including in layers that 
do not receive extrastriate feedback and do not have long- 
range horizontal connections. Both phenomena have been ob- 
served in the lateral geniculate nucleus (LGN) , and are likely 
to be partially inherited by VI cells via feedforward connec- 
tions from LGN 9,16 . Finally, there is experimental evidence 
for contextual modulations mediated by local short-range con- 
nections in cat . 

In this paper we show that local short-range connections 
in VI and LGN input are, in principle, sufficient to explain 
these two extraclassical receptive field phenomena in layers 
AC a and AC (3 of macaque. We do this by means of a large- 
scale computational model which is constructed, as much as 
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possible, from basic, established experimental data. We sug- 
gest neural mechanisms for the phenomena by analyzing the 
synaptic inputs that generate them. An illustration of the 
model's architecture is given in Figure 1. A brief summary of 
the model is given in Methods. 

Methods 
The Model 

Our model consists of 8 ocular dominance columns and 64 
orientation hypercolumns (i.e. pinwhccls), representing a 16 
mm 2 area of a macaque VI input layer 4Ca or 4C/3. The 
model contains approximately 65,000 cortical cells and the 
corresponding appropriate number of LGN cells. Our cor- 
tical cells are modeled as conductance based integrate-and- 
fire point neurons, 75% are excitatory cells and 25% are in- 
hibitory cells. Our LGN cells are rectified spatio-temporal 
linear filters. The model is constructed with isotropic short- 
range cortical connections (< 500/xm), realistic LGN receptive 
field sizes and densities, realistic sizes of LGN axons in VI, 
and cortical magnification factors and receptive field scatter 
that are in agreement with experimental observations. We 
will only give a very brief description of the model here, it is 
explained in detail in Supplementary Materials. Some back- 
ground information can also be found in previous work 18 ' 19 
by one of the authors (JW). 

Dynamic variables of a cortical model-cell i are its mem- 
brane potential Vi(t) and its spike train Si(t) — ^ fc 5(t — ti,k), 
where t is time and t^k is its fcth spike time. Membrane po- 
tential and spike train of each cell obey a set of N equations 
of the form 

dv ' 

C i~^ = -9L,i{Vi-V L ) - gE,i(t,[S] E ,T]E){Vi -v E ) 

[S]i,Vi){vi - vi) , i = 1, . . . , N . (1) 

These equations are integrated numerically using a second 
order Runge-Kutta method with time step 0.1 ms. Whenever 
the membrane potential reaches a fixed threshold level vt it 
is reset to a fixed reset level vr and a spike is registered. The 
equation can be rescaled so that Vi(t) is dimensionless and 
d = 1, v L = 0, v E = 14/3, vi = -2/3, v T = 1, v R = 0, and 
conductances (and currents) have dimension of inverse time. 

The quantities g E ,i(t, [S],t]e) and gi y i{t, [S],r]i) are the ex- 
citatory and inhibitory conductances of neuron i. They are 
defined by interactions with the other cells in the network, 
external noise t]e(i), and, in the case of g E ,i possibly by LGN 
input. The notation [«S]b(i) stands for the spike trains of all 
excitatory (inhibitory) cells connected to cell i. Both, the 
excitatory and inhibitory populations consist of two subpop- 
ulations Vk(E) and Vk{I), k = 0,1, a population that re- 
ceives LGN input (k = 1) and one that does not (k = 0). 
In the model presented here 30% of both the excitatory and 
inhibitory cell populations receive LGN input. We assume 



noise, cortical interactions and LGN input act additively in 
contributing to the total conductance of a cell, 

g E ,i(t, [S]e,ve) = mM + gZ^ We) + 5 l9 f GN {t) 

gili (t, [5]!,*?/) = ViAt)+9iJ(t, [S]!) , (2) 

where <5; = I for i e {V e (E),V e {I)}, £ = 0,1. The terms 
g??J{t, [<S] M ) are the contributions from the cortical excitatory 
(/i = E) and inhibitory (/j = I) neurons and include only 
isotropic connections, 



/ + °° ds E E C k ^{\\x i -x j \\)G,At-s)S j {s), (3) 
J -°° k=o jev k (n) 

where i € Vk'(fjf) Here Xi is the spatial position (in cortex) of 
neuron i, the functions G MJ (r) describe the synaptic dynam- 
ics of cortical synapses and the functions C ,' (r) describe the 
cortical spatial couplings (cortical connections). The length 
scale or excitatory and inhibitory connections is about 200^m 
and lOO^tm respectively. 

An important class of parameters are the geometric param- 
eters, which define and relate the model's geometry in visual 
space and cortical space. Geometric properties are different 
for the two input layers 4Ca,(3 and are different at different 
eccentricities. As said, the two extraclassical phenomena we 
seek to explain are observed to be insensitive to those differ- 
ences 4 ~ 6 ' 8 . In order to verify that our explanations are con- 
sistent with this observation, we have performed numerical 
simulations for four different sets of parameters, correspond- 
ing to the 4Ca, f3 layers at para-foveal eccentricities < 5° and 
at eccentricities around 10°. These different model configu- 
rations are referred to as M0, M10, and P0, P10 in the text. 
Reported results are qualitatively similar for all four configu- 
rations unless otherwise noted. 

In agreement with experimental findings (see references in 
18 ), the LGN neurons are modeled as rectified center-surround 
linear spatiotemporal filters. A cortical cell, j € V\{^) is 
connected to a set N^ N of left eye LGN cells, or to a set 

N k G j N of ri g ht e y e LGN cclls > 

gf GN (t)= ]T [g° l+ gY 

r ds [ d 2 y Gp GN (t - 8 )C t { \\y e - y\\) I(y, *)]+ . (4) 

J — oo J 

where Q = L or R. Here [a;]+ = x if x > and [a;]+ = 
if a; < 0, Ci(r) and Gf GN (r) are the spatial and temporal 
LGN kernels respectively, ye is the receptive field center of the 
£th left or right eye LGN cell, which is connected to the jth 
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cortical cell, I(y, s) is the visual stimulus. The parameters 
g\ represent the maintained activity of LGN cells and the 
parameters gj measure their responsiveness to visual stimuli. 
Their numerical values are taken to be identical for all LGN 
cells in the model, g® = 2s _1 and gj = 25 cd~ 1 m 2 s~ 2 . The 
LGN kernels are of the form 20 
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and 



Ct (r) = ±(1 - K t 
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,-t/t 2 ) 



if T < T° 
if T > t" 



(5) 



-(r/o 



-(r/a s , e f 



(6) 

where k is a normalization constant, <j c j and a Sl e are the 
center and surround sizes respectively, and Kt is the in- 
tegrated surround- center sensitivity. The temporal kernels 
are normalized in Fourier space, \Gf GN (uj)\du = 1, 

Gf GN {cj) = (27T)- 1 f^Gf^i^e-^dt. For the magno 
cases (MO, M10) the time constants t\ = 2.5 ms and T2 = 7.5 
ms and c = (ti/t 2 ) 6 so that Gf GN (Q) = 0, in agreement 
with experiment 20 . For the parvo cases (P0, P10) the time 
constants t\ = 8 ms and T2 = 9 ms and c = 0.7(ti/t2) . 
The delay times t° are taken from a uniform distribution 
between 20 ms and 30 ms, for all cases. Sizes for center 
and surround were taken from experimental data 21-25 and 
were a c ^ t = a c = 0.1°, 0.2°, 0.04°, 0.0875° (centers) and 
a s , e = ct s = 0.72°, 1.4°, 0.32°, 0.7° (surrounds), for M0, M10, 
P0 and P10 respectively. The integrated surround-center sen- 
sitivity was in all cases Kg = 0.55 25 . By design, no diversity 
has been introduced in the center and surround sizes in order 
to demonstrate the level of diversity resulting purely from the 
cortical interactions and the connection specificity between 
LGN cells and cortical cells (i.e. the sets Nq GN , see speci- 
fications below). Further, no distinction was made between 
ON-center and OFF-center LGN cells other than the sign re- 
versal of their receptive fields (± sign in Eq. |HJ). The LGN 
RF centers yi were organized on a square lattice with lat- 
tice constants cr c /2, a c , <r c /2, and a c /2 for M0, M10, P0 and 
P10 respectively. These lattice spacings and consequent LGN 
receptive field densities imply LGN cellular magnification fac- 
tors that are in the range of the experimental data available 
for macaque 26 ' 27 . The connection structure between LGN 
cells and cortical cells, given by the sets Nq GN , is made so 
as to establish ocular dominance bands and a slight orienta- 
tion preference which is organized in pinwheels 28 . It is further 
constructed under the constraint that the LGN axonal arbor 
sizes in VI do not exceed the anatomically established val- 
ues of 1.2 mm for magno and 0.6 mm for parvo cells 29 ' 30 . A 
sketch of the model is given in Figure 1. Further details are 
given in Supplementary Information. 

Some of the geometric differences (in the model as well 
as in the true biological situation) may be expressed by the 



dimensionless parameter Q — 



where 
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cortical magnification factor, a c is the LGN receptive field 
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Figure 1: Model architecture, (left top) A typical cluster of ON 
(blue circles) and OFF (red dots) M-LGN cells that feed into one 
cortical cell. Receptive field centers of LGN cells are organized 
on a square lattice (orange), (left bottom) Some typical M-LGN 
axons in our model- VI. Points of the same color are cortical cells 
that connect to the same LGN axon, (right) Pinwheel structure 
and ocular dominance columns for M10 model, constructed from 
averaged responses in the spirit of optical imaging experiments 28 . 
All sample cells used to study extraclassical phenomena are taken 
from within the white dashed rectangle (see Supplementary Infor- 
mation for details). 



size (center size), and £ c is a characteristic length scale for 
the excitatory cortical connectivity. Substituting numerical 
values taken from experimental data, this parameter is 1, 0.57, 
0.4, and 0.25, for M0, M10, P0 and P10 respectively. At 30° 
eccentricity, the experimental data suggests values for this 
parameter not very different from its values at 10° (f2 = 0.5 
for M30 and fi = 0.25 for P30). 

Visual stimuli and data collection 

The stimulus used in this paper to analyze the phenomena of 
surround suppression and receptive field growth at low con- 
trast is a drifting grating confined to a circular aperture, sur- 
rounded by a blank (mean luminance) background. The lu- 
minance of the stimulus is given by I{y,t) = Ia(l + ecos( ujt — 
k-y + (f>)) for < r A and I(y,t) = I for \\y\\ > r A , with av- 
erage luminance Iq, contrast e, temporal frequency u>, spatial 
wave vector k, phase <f>, and aperture radius ta- The aper- 
ture is centered on the receptive field of the cell and varied 
in size, while the other parameters are kept fixed and set to 
preferred values. All stimuli are presented monocularly. As 
the aperture size increases the response of a VI cell to such 
stimuli typically reaches a maximum, after which it settles 
down to some steady level. The aperture size for which the 
response reaches its maximum is sometimes referred to as the 
"classical" receptive field size 4 > 10 - 31 . We will simply refer to 
the minimum aperture radius for which the response 
is > 95% of its maximum as the receptive field size (r). We 
define the surround size (R) as the minimum aperture radius 
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> r for which the suppression / s (r^) = fmax — f( r A) is > 95% 
of its maximum. We define the asymptotic response /oo as the 
average response beyond R. We define the suppression index 
(SIi) as the relative surround suppression, 



Sh = 



fn 



fmax fo 



(7) 



where fo is the response to a blank stimulus. The suppression 
index SIi is similar to the one used in 8 , but different from 
the integrated suppression index used in 4 . 

The primary data, i.e. responses and conductances as a 
function of aperture size for single eye stimulation, are ob- 
tained for samples of approximately 200 cells for each config- 
uration, containing about an equal number of simple and com- 
plex cells. Each stimulus was presented for 3 s and preceded 
by a 1 s blank stimulus. The procedure was repeated five 
times with different initial conditions and noise realizations. 
Standard errors in cycle-trial average responses and conduc- 
tances are negligibly small. The experiments were performed 
at "high" contrast, e = 1, and "low" contrast, e = 0.3. More 
precise definitions and further details are in Supplementary 
Information. 

DOG & ROG models 

In the Difference-of-Gaussians (DOG) model 4 ' 6 ' 10 , the re- 
sponse /(ta) is fit to 

f(r A ) = 



\/7T 



dy - Kj 



dy 



(8) 



In this model the response is assumed to arise from a summa- 
tion of background activity (fo), "excitation" (spatial scale 
<je) and "inhibition" (spatial scale 07). The integrated sup- 
pression index SI2 is defined as 

Kicri 



Sh 



Ke<je 



(9) 



As is true for SIi, 5/2 can be larger than one, indicating 
surround suppression beyond the background response. 

Given the validity of the rectification model (Eq. I13fl . we 
can "derive" the DOG model by the substitutions qe,i ~ 
J Q rA e~{ y l aE ^ dy. Obviously, identification of the terms in 
the DOG model with the actual excitatory and inhibitory 
inputs can be little more than symbolic. 

The Ratio-of-Gaussians (ROG) model 8 is defined by 



f(r A ) =fo + 



2 rr* 

sjn Jo 



e -(y/wcf dy y 



1 



(_2_ f 



c 



2 dy 



(10) 



In this model the response beyond the background response 
is assumed to arise from a division of center activity ("exci- 
tation", spatial scale w c , gain k c ) and surround activity ("in- 
hibition", spatial scale w s , gain k s ). 



From the slaving of the membrane potential (Eq. I12fl . it is 
a simple matter to "derive" the ROG model from the stan- 
dard rectification model. Equation (|12[1 can be rewritten such 
that the numerator (N) and denominator (D) represent a rec- 
tified weighted difference of the excitatory and inhibitory con- 
ductances, and the total conductance gr, respectively. The 
ROG model used in 8 is then obtained by the substitutions 



N,D-1 



TA - 

e 



dy 



. As is also true here, identifi- 



cation of the terms in the ROG model with actual excitatory 
and inhibitory inputs, can be little more than symbolic. 

Results 

Classical response properties 

Classical responses do not specifically address size effects of 
the stimulus. These are response properties such as orien- 
tation tuning, spatial and temporal frequency tuning, distri- 
bution of response modulations for drifting grating stimuli 
(simple & complex cells) etc.. One of our model's strong 
accomplishments is that it produces, with the same fixed pa- 
rameters, a wide range of classical responses as well as the 
two extraclassical responses which are the focus of this paper. 

Classical response properties are important because they 
set the context in which extraclassical responses occur. First, 
because extraclassical responses (responses evoked from out- 
side the classical receptive field) are not know to occur with- 
out sufficient stimulation of the classical receptive field. Sec- 
ond, extraclassical responses, in particular, depend strongly 
on how the cell's environment in the cortex is responding. 
The responses of the cells that make up this environment will 
display an enormous diversity to any particular fixed stimu- 
lus. A cell's relevant cortical environment is generally made 
up of cells with, for instance vastly different orientation, spa- 
tial and temporal tuning widths and preferences. A reason- 
able response of a cell's environment is thus accomplished if 
a model's classical responses are realistic, i.e. agree with ex- 
perimental data. 

Some classical response properties of the model are illus- 
trated in Figures 2 & 3. All plots are for the MO configura- 
tion (see Methods) but the other configurations yield similar 
results. The sharpness and spatial distribution of orientation 
tuning in the model is illustrated in Figure 2A. As a measure 
for orientation tuning we used the circular variance (CV) , 



CV = 1- 



/ r(6) exp(2i8)d6 



Jr(9)d6 



(11) 



Here r(0) is the mean firing rate and 8 the grating angle. The 
smaller the CV, the sharper the orientation tuning. Cells with 
CV = respond at just one angle, and hence are very sharply 
tuned. Cells with CV — 1 respond identically at all angles, 
and hence are not tuned for orientation. In Figure 2A we color 
coded the CV for all cells within the white dashed rectangle 
in Figure 1 of the paper. The stimulus (drifting grating) was 
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presented to one eye, the other eye received no visual input 
(sec Materials and Methods for details). Pixels colored black 
indicate cells that do not show a significant response under 
visual stimulation and are mostly cells that receive their input 
from the other eye. Notice that our model cortex is filled with 
sharply tuned cells, moderately tuned cells, and untuned cells, 
as is the primary visual cortex of macaque. Notice also that 
there is no particular spatial organization of the sharpness of 
orientation tuning. 
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Figure 2: Some classical response properties (selectivity and tun- 
ing) of the M0 version of the model, other configurations yield 
similar results, (left) Spatial distribution of the sharpness of orien- 
tation tuning, expressed in the circular variance (CV, color coded), 
for cells within the white rectangle in Figure 1 of paper. Black pix- 
els are cells that do not show sufficient response under single eye 
stimulation, (right) Temporal and spatial frequency tuning curves 
of some M0 sample cells and histograms of preferred temporal and 
spatial frequencies for cells in the M0 sample set. Thick black 
curves refer to the LGN cells. 

Our model yields a realistic diversity in spatial and tempo- 
ral frequency tuning properties. This is illustrated in Figure 
2B-E. In observing this diversity in the model it is impor- 
tant to note that all LGN cells in a particular configuration 
are identical by construction. Their spatio-temporal tuning 
properties are indicated by the thick black curves in (B) and 
(D). Notice that, as in reality, the preferred temporal fre- 
quencies of our cortical cells are smaller than the preferred 
temporal frequency of our LGN cells. The main reason for 
this is the inclusion of slow (NMDA) excitation in the model. 
The diversity seen in the spatial frequency tuning of our cor- 
tical cells, particularly in the bandwidth, is mostly cortical in 
origin but also partially results from the spatial diversity in 
the feedforward (LGN) connections. 

Our first example of classical spatial summation properties 
of the model is provided in Figure 3A, B. Shown are aver- 
aged response waveforms of spike train and membrane po- 
tential in response to a drifting grating in Figure 3A. These 
are responses of a simple and a complex cell in the model 
for several angles of the grating at their preferred spatial and 
temporal frequencies. These modulations in the spike train 
at the preferred angle are frequently used to classify simple 



and complex cells in VI. A cell is called complex whenever 
F1/F0 < 1, and simple otherwise, where FI is the first har- 
monic of the response and FO the average. The distribution 
of the spike train modulation index F1/F0 over our cell pop- 
ulation is shown in Figure 3B (top left). Our model cortex 
contains about an equal amount of simple and complex cells. 
The bimodal shape of the distribution of the spike train mod- 
ulation index agrees well with experimental data 32 . 
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Figure 3: Classical spatial summation properties of the P0 ver- 
sion of the model, other configurations yield similar results. (A) 
Response waveforms for a simple and a complex cell in response to 
a drifting grating, for a number of different angles. (B) Distribu- 
tions (normalized to peak value one) of modulation index (F1/F0) 
in spike train (top left) and membrane potential (bottom left), and 
spatial distribution of the modulation index (spike train) for cells 
within the white rectangle in Figure 1 of the paper. Black pixels are 
cells that do not show sufficient response under single eye stimula- 
tion (right). (C) Response waveforms for a simple and a complex 
cell in response to a contrast reversal grating at the preferred angle, 
for a number of different spatial phases. (D) Distributions of the 
phase averaged F2/F1 ratio for spike train (top) and membrane 
potential (bottom) for responses to a contrast reversal grating at 
the preferred angle. 

It is easy to understand how the diversity in response modu- 
lations occurs in our model. The modulations enter our model 
cortex via the LGN input received by 30% of the cortical cells. 
The phases of these LGN inputs into the different cortical cells 
vary randomly on [0, 27r]. This is so because of the receptive 
field off-sets of the clusters of LGN cells connected to different 
cortical cells, the difference in shape (symmetry) of the clus- 
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ters themselves, and the diversity in temporal delays in the 
LGN kernels. A cell receives input from many other cells, thus 
a cell's excitatory and inhibitory inputs will show stronger or 
weaker modulations dependent on its specific environment in 
the network and whether or not it receives LGN input. Inter- 
play between the strengths and phases of the modulations in 
these inputs and cell specific parameters ultimately determine 
the modulation in the cell's spike and membrane potential re- 
sponse. Most cells that receive LGN input are simple cells 
(80% in our model) and most cells that do not receive LGN 
input are complex (70% in our model). As mentioned earlier, 
strength parameters have been set so that the distribution of 
these modulations in the spike responses is in agreement with 
experimental data. 

The distribution of modulations in the membrane poten- 
tial (with respect to a blank stimulus) is shown in Figure 3B 
(bottom left). Notice that the bimodality present for spike 
train modulations is not present in the modulation index dis- 
tribution for the membrane potential. However, our model 
predicts (not shown explicitly) that the classification of sim- 
ple and complex cells can equally well be made in terms of 
the membrane potential modulation index distribution, the 
two modes in this case being its "core" (complex cells) and 
its "tails" (simple cells). Also notice that our model predicts 
a "gap" in the distribution at small negative values. The dis- 
tribution of modulations in the membrane potential has not 
yet been observed experimentally in macaque. Some data for 
cat has recently been published 33 and they do not contradict 
the predictions based on our model. The spatial distribution 
of the modulation index (spike train) across all cells within 
the white dashed rectangle in Figure 1 of the paper is shown 
in Figure 3B (right). The stimulus was again presented to one 
eye, and pixels colored black indicate cells that do not show 
a significant response under visual stimulation. The figure 
shows that simple and complex cells are randomly distributed 
across space in our model cortex, i.e. there is no particular 
spatial organization of the modulation index. 

Another example of classical spatial summation properties 
in our model is provided in the remainder of Figure 3. Aver- 
aged response waveforms of spike train and membrane poten- 
tial in response to a contrast reversal grating at the preferred 
angle are shown in Figure 3C. Shown are the responses of 
a simple and a complex cell in the model for several spatial 
phases of the grating. Simple cells perform an approximately 
linear spatial summation, that is, their responses contain a 
dominant Fl component and the spatial phase dependence 
of their response waveform is similar to the spatial phase 
dependence of the stimulus. Complex cells respond nonlin- 
early, their response waveform is relatively insensitive to spa- 
tial phase and contain a dominant (frequency doubling) F2 
component. The distribution of the phase-averaged F2/F1 
for our model cells is shown in Figure 3D. For what concerns 
the spike train waveforms (top), this quantity displays a weak 
bimodality and its behavior for our model cells agrees with 
experimental data 34 , complex cells having mostly F2/F1 > I 



and simple cells F2/F1 < 1. Note that this property of our 
model cells follows naturally, without any parameter adjust- 
ments, after the strength parameters have been set to achieve 
essentially only orientation tuning and a proper distribution 
of response modulations in response to a drifting grating (Fig. 
3B), as mentioned earlier (see also 19 ). 

It is easy to understand how the diversity in F2/FI occurs. 
As explained in 19 , for a contrast reversal grating stimulus 
each total LGN input into a cortical cell has in general a 
dominant Fl component with a phase close to either or n, 
determined by the relative positions of the ON and OFF sub- 
fields of the corresponding cluster. The cortical excitatory 
and inhibitory inputs in a cell will thus have a relative strong 
F2 component since they arise from many other cells. The 
actual strengths of Fl and F2 components in a cell's excita- 
tory and inhibitory inputs thus depends on the cell's specific 
environment in the network and on whether it receives LGN 
input or not. Interplay of these inputs and cell specific pa- 
rameters determine the F2/F1 ratio in the cell's spike and 
membrane potential response. Clearly, most cells that receive 
LGN input (simple) will have F2/F1 < 1 and most cells that 
do not receive LGN input (complex) will have F2/F1 > 1. 

No experimental data is available for the distribution of 
F2/F1 of the membrane potential waveforms. Our model's 
prediction is shown in Figure 3D, bottom. Our model pre- 
dicts that, quite contrary to the situation for the modula- 
tion index F1/F0, the (weak) bimodality of the distribution 
of F2/F1 for spike waveforms is not eliminated, but rather, 
becomes more pronounced in the F2/F1 distribution for mem- 
brane potential waveforms. This, in fact, can be understood 
quite simply from a standard rectification model, in which 
the membrane potential waveforms are subjected to a thresh- 
old to give the spike waveforms. For complex cells, both the 
membrane potential and spike responses will contain a strong 
F2 component. Hence in this case practically all of the mem- 
brane potential waveform will be above threshold, so that 
evaluation of F2/F1 will yield about the same result for spike 
waveforms as for membrane potential waveforms. This is also 
apparent in Figure 3D: the F2/F1 > 1 section of the two dis- 
tributions (in top and bottom panels) is similar. For simple 
cells, the membrane potential and spike responses will contain 
a dominant Fl component, and for both responses about an 
equally small F2 component. Because of the rectification, the 
Fl component in the membrane potential waveform is sub- 
stantially reduced in the spike waveform. Hence, F2/F1 will 
turn out substantially smaller when evaluated for membrane 
potential responses. This is again apparent in Figure 3D: the 
F2/F1 < 1 (simple cells) section of the distribution for spike 
trains (top) is shifted to the left in the F2/F1 distribution for 
membrane potential waveforms (bottom). 

Extraclassical spatial summation 

In this section we summarize the extraclassical results for our 
model and compare them with experimental data. Quite con- 
trary to classical response properties, the two extraclassical 
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responses we discuss in this paper are, as in experimental 
data, not substantially different for simple and complex cells 
and results that follow are not type specific. An example of 
the surround suppression and contrast dependent receptive 
field size observed in our model is shown in Figure 4A, for 
a cell from the para-foveal AC a model (MO) (see Methods). 
Shown are responses for both firing rate and membrane po- 
tential, at high (solid) and low (dashed) contrast. 

Distributions of receptive field and surround sizes for the 
4C/3, 10°-eccentricity model (P10) are shown in Figure 4B 
& C. The distributions for the other model configurations are 
given in Supplementary Information. Receptive field sizes and 
surround sizes in our model show excellent agreement with 
experimental data 6,8 . This is true for the mean values, for 
the diversity, as well as for their dependence on eccentricity 
from para-foveal to 10° and 30° eccentricity 8 . 



Mechanisms of surround suppression 

The DOG and ROG models arc phenomenological models and 
do not provide much insight into the neural mechanisms of the 
phenomena. Both models miss an essential feature of the ex- 
citatory and inhibitory inputs, which is that these inputs gen- 
erally themselves show surround suppression 7 . In our model 
we similarly observe a significant suppression in both conduc- 
tances, an example is shown in Figure 5C. This cell shows 
that, unlike suggested by the DOG and ROG picture, sur- 
round suppression in the spike response takes place entirely in 
the region of decreasing synaptic inputs (conductances). We 
can say that the surround suppression of this cell is caused 
by a decrease of excitation, since the decrease of inhibition 
could not by itself suppress the cell's response. This cell is 
not atypical in our model and the above scenario is indeed 
how surround suppression works in about 50% of the cells. 



The distribution of surround suppression and receptive field 
growth for the MO configuration of our model is given in Fig- 
ure 4D & E. In agreement with experimental data, the shape 
of the distribution of the suppression index SI\ (Methods) 
is skewed to low suppression 8 . (Cells without surround sup- 
pression have SI\ = 0, cells with fully suppressed response for 
large stimuli have SIi = 1.) We also see a slight increase of 
the average suppression for low contrast 4,8 . The average sup- 
pression index (all eccentricities, see Methods) is SIi ~ 0.2 
and this is about half of what is observed experimentally 8 . 
The receptive field and surround growths (Fig. 4E) are ex- 
pressed as ratios, r_/r + and R_/R + respectively. We observe 
an average growth by about a factor of two in both receptive 
field size (r_/r + ~ 2) and surround size (R-/R+ ~ 2). This 
receptive field growth is a little less than what is observed in 
experiments 5,s . 

We fitted our data with the Difference-of-Gaussians (DOG) 
and Ratio-of-Gaussians (ROG) 4,8 models (Supplementary In- 
formation). We obtain for the integrated suppression index 
SI2 ~ 0.4. Average growth ratios for the excitatory space 
constant is cr^/er^ ~ 1.5 (both DOG and ROG, all eccentric- 
ities). The suppression index and growth ratios are again less 
than what is seen experimentally (0.6 and 2.3 respectively 4 ). 
In agreement with experimental studies 4, 8 there is no system- 
atic dependence of suppression on contrast in either index SIi 
or 5/2. This observation is illustrated in Figure 5 A. 

All of the above findings are based on spike responses. 
Membrane potential responses yield qualitatively similar re- 
sults, but, due to the spike threshold, suppression in the mem- 
brane potential is systematically smaller. This is illustrated 
in Figure 5B. The same observation has also been made ex- 
perimentally in cat 7 . 

A more extensive summary of our model data for different 
eccentricities and including receptive field sizes and surround 
sizes is given in Supplementary Information. 
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Figure 4: Summary of extraclassical spatial summation in the 
model. (A) Response as function of aperture sized for a cell from 
the M0 model (see Methods). Shown are firing rate (black) and 
membrane potential (gray) for high (squares) and low (circles) con- 
trast. Standard errors are negligibly small. (B & C) Receptive 
field and surround sizes for the P10 model at high (unfilled) and 
low (shaded) contrast. The diversity of responses produced by the 
model is similar to what is seen in experimental data 6,8 . (D) Dis- 
tribution for the M0 model of the suppression index 5*11 at high 
(unfilled) and low (green shaded) contrast. All suppression is ex- 
clusively due to short-range cortical connectivity. (E) Distributions 
for the M0 model of the ratios of the receptive field and surround 
sizes at low and high contrast, r-_/r+ (blue shaded) and R-/R+ 
(unfilled). (Wilcoxon test on ratio larger than unity: p < 0.001 for 
both receptive field and surround growth). For a more complete 
summary of our model data see Supplementary Information. 
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Analysis of the surround suppression in our model is based 
on the fact that the average membrane potential (i>fc(i, r A )) 
and instantaneous firing rate (S k (t, r A )) (of the fc-th neuron) 
are well-approximated by 

(v k (t,r A )) « V k = (12) 
K9T.k) 

and 

(S k (t,r A )) « f k = 6 k {(I D , k )-(g T ,k}-Ak} + ■ (13) 

Here [x] + = a; if x > and [x]+ = if x < 0, and for most 
cells good approximations are obtained with a gain S k and 
threshold A k that do not depend on the aperture radius r A 
nor time. The total conductance gT,k and difference current 
In,k are given by 

9r,k = 9l+ 9E,k + 9i,k (14) 

lD,k = 9E,k V E - gi, k \Vi\ . (15) 

Equations (|12[1 and {T5)| allow us to base our analysis directly 
on the (cycle-trial averaged) conductances as a function of 
the aperture radius r A and time. In what follows we drop 
the averaging notation (•), assuming it unless stated other- 
wise. Given Equations (|12(l and Ijl3(l there are three ways that 
surround suppression of spike train and membrane potential 
could arise, namely (A) dgE.k/dr A > and dgx,k/dr A > 0, 
(B) dg E ,k/dr A < and d 9l , k /dr A < 0, or (C) dg E , k /dr A < 
and dgi, k /dr A > 0. In other words, surround suppression is 
caused by (A) an increase in the inhibitory conductance, or 
(B) a decrease in the excitatory conductance, or (C) both (A) 
and (B) simultaneously. 

Examples of this analysis for a (simple) cell receiving LGN 
input and a (complex) cell that does not receive LGN in- 
put are given in Figure 6. The cycle-trial averaged conduc- 
tances for apertures around the aperture of maximum re- 
sponse (marked by an asterisk) are shown in Figure 6C-F. For 
example, by comparing the conductances for aperture "aster- 
isk" and the aperture for which the suppression is completed, 
we see that at high contrast the suppression mechanism for 
the simple cell is (A) and for the complex cell it is (B). At 
low contrast the suppression mechanisms are (C) and (B) re- 
spectively. 

We observe all three mechanisms A, B, and C in our model. 
Typically, different mechanisms act sequentially as the aper- 
ture size r A increases from receptive field size r to surround 
size R, while in some cases we find that different mechanisms 
are active during different times in the stimulus cycle. 

As may be clear from Figure 6, identifying the mechanisms 
for surround suppression based on Eq. H12|) and l|13l) can be 
rather more subtle than just comparing the mean (F0) con- 
ductance, its first harmonic (Fl) or the peak conductance 
(~F0+F1). However, we find that for most cells, an analy- 
sis using the F0+F1 components of the conductances allows 
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Figure 5: Two example cells, an MO simple cell which receives 
LGN input (left) and an M10 complex cell which does not receive 
LGN input (right). (A & B) Responses as function of aperture size. 
Mean responses are plotted for the complex cell, first harmonic for 
the simple cell. Apertures of maximum of responses (i.e. recep- 
tive field sizes) are indicated with asterisks (red=high contrast, 
green=low contrast). (C & D) Conductances for high contrast at 
apertures near the maximum responses. (E & F) Conductances 
for low contrast at apertures near the maximum responses. Panels 
C-F each consist of seven sub-panels giving the cycle-trial averaged 
conductances as function of time (relative to cycle) and aperture 
size. Asterisks indicate corresponding apertures of maximum re- 
sponse in A-B. 

identification of the suppression mechanisms. Comparing con- 
ductances at r A = r and at r A = R in this way, we find that 
at low contrast all three mechanisms are about equally preva- 
lent, while at high contrast mechanism A is somewhat more 
likely than B and C. 

Mechanisms of contrast dependent receptive 
field size 

The DOG model suggests that growth in receptive field size 
at low contrast is due to an increase of the spatial summation 
extent of excitation 4 (<te). This was partially confirmed ex- 
perimentally in cat primary visual cortex 7 . Although it has 
been claimed 8 that the ROG model would explain receptive 
field growth solely from a change in the relative gain param- 
eter k s , we believe this is incorrect. Since there is a one-one 
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Figure 6: Relations between some key response measures for the 
MO configuration of the model, other cases yield qualitatively sim- 
ilar results. (A) Scatterplot of surround suppression at low and 
high contrast expressed in the two different suppression indexes 
5/1 (black) and SI2 (orange). (B) Scatterplot of surround sup- 
pression in spike train and membrane potential at high (black) 
and low (gray) contrast. (C) Spike responses, membrane potential 
responses, and cortical conductances as function of aperture size 
for a model cell which shows about 50% surround suppression (in 
spike train). Notice the surround suppression of the conductances. 



relation between k s and the surround suppression, this would 
imply that contrast dependent receptive field size simply re- 
sults from contrast dependent surround suppression, which 
contradicts experimental data 4,8 . Thus, as does the DOG 
model, the ROG model predicts that contrast dependent re- 
ceptive field size is due to contrast dependence of the spatial 
summation extent of excitation. As we show below, our simu- 
lations confirm an average growth of spatial summation extent 
of excitation (and inhibition) at low contrast. However, this 
growth is neither sufficient nor necessary to explain receptive 
field growth. 

From Eq. i|12[) and (|13[) it follows that a change in receptive 
field size in general results from a change in behavior of the 
relative gain parameter, defined as 



G(r A ) = 



dg E /dr A 
dgijdr A 



Note that this is a rather different parameter than the "sur- 
round gain" parameter k s used in the ROG model. (For ex- 
ample, unlike k s , G(r A ) is not simply related to the degree of 
surround suppression.) Qualitatively, the conductances show 
a similar dependence on aperture size as the membrane po- 
tential responses and spike responses in that they display sur- 
round suppression (Fig. 5C). Receptive field sizes based on 
these conductances are a measure of the spatial summation 
extent of excitation and inhibition. 

A change in the spatial summation extent of §e and/or gi 
is just one of the many ways to change the behavior of G and 
consequently the receptive field size. For example, some other 
possibilities are illustrated by the two cells in Fig. 6. These 
cells show, both in spike and membrane potential responses, 
a receptive field growth of a factor of 2 (left) and 3 (right) at 
low contrast. However, for the left cell, the spatial summation 
extent of excitation at low contrast is one aperture less than 
at high contrast, and for inhibition at low contrast it is one 
aperture larger than at high contrast. For the cell in the 



right the spatial summation extents of both excitation and 
inhibition do not change with contrast. 

In a similar way as for spike train responses, we also ob- 
tained receptive field sizes for the conductances. As shown in 
Figure 7, both excitation and inhibition also show on the aver- 
age an increase in their spatial summation extent as contrast 
is decreased, but the increase is in general smaller than what 
is seen for spike responses, particularly for cells that show sig- 
nificant receptive field growth, say (r+/r_ > 1.5). (Wilcoxon 
test on ratio of growth ratios larger than unity: p < 0.05 
(all cells, excitation, Fig. 7B), p < 0.15 (all cells, inhibition, 
Fig. 7C), p < 0.001 (cells with receptive field growth rate 
7' + /r_ > 1.5, both excitation and inhibition.) Further, Fig- 
ure 7B & C show that, although some increase in the spatial 
summation extent of excitation and inhibition is in general the 
rule, this increase is rather arbitrary and bears not much rela- 
tion with the receptive field growth based on spike responses. 
For cells in the sample with larger receptive field growths (fac- 
tor of ~ 2 and greater) this growth is always considerably less 
than the growth based on their spike responses. The same 
conclusions follow from membrane potential responses (not 
shown) . 

A more precise analysis based on the relative gain param- 
eter is given in Supplementary Information. For cells with 
significant receptive field growth, (r + /r_ > 1.5) we are able 
to identify a systematic mechanistic property. We find that 
for more than 50% of such cells, a transition takes place from 
a high contrast RF size less or equal to the spatial summation 
extent of excitation and inhibition, to a low contrast receptive 
field size which exceeds both. 
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Figure 7: (A) Joint distribution of high and low contrast recep- 
tive field sizes, r+ and r_, based on spike responses. All scales 
are logarithmic, base 10. All distributions are normalized to peak 
value one. Receptive field growth at low contrast is clear. Av- 
erage growth ratio is 1.9 and is significantly greater than unity 
(Wilcoxon test, p < 0.001). (B & C) Joint distributions of re- 
ceptive field growth and growth of spatial summation extent of 
excitation (B) and inhibition (C) (computed as ratios). There is 
no simple relation between receptive field growth and the growth 
of the spatial summation extent of excitatory or inhibitory inputs. 
For cells in the sample with larger receptive field growths (factor 
of ~ 2 and greater) this growth is always considerably larger than 
the growths of their excitatory and inhibitory inputs. 
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LGN contributions 

The extraclassical responses in our model discussed so far are, 
by construction, exclusively the result of cortical interactions 
and not inherited from LGN inputs. This is due to our use of 
the standard center-surround model for LGN receptive fields 
(Supplementary Information). Thus our LGN cells show nei- 
ther surround suppression nor contrast dependent receptive 
field size. For what concerns the latter, this is true at all 
spatial frequency. Surround suppression, however, is absent 
at the relevant (optimal) spatial frequencies for our cortical 
cells (SIi(gig n ) < 0.01), but does appear at lower spatial fre- 
quencies. We used this fact to study the transfer of LGN 
surround suppression to cortical cells. Surround suppression 
of our LGN cells for different spatial frequencies is illustrated 
in Figure 8A. No suppression occurs for cortical optimal spa- 
tial frequencies kc, suppression starts at spatial frequencies 
of about 0.5fcc, it becomes stronger for smaller spatial fre- 
quencies and is about 25% at 0.25kc- Further, we find that 
at 0.25fcc the surround is not yet able to evoke responses on 
its own, and the suppression thus in some sense qualifies as 
"extraclassical" . 

(A) (B) (C) 
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Figure 8: Transfer of LGN surround suppression to cortical cells. 
(A) Responses (Fl) of the P10 model LGN cells as function of aper- 
ture size, at different spatial frequencies (all LGN cells in a par- 
ticular configuration are identical). Average cortically preferred 
spatial frequency is kc (red), at which our LGN cells show no 
surround suppression. The LGN cells show about 25% (classical) 
surround suppression at a spatial frequency k — 0.25fcc (blue). (B 
& C) Surround suppression and prevalence of suppression mech- 
anisms at spatial frequencies k = kc (unfilled) and k — 0.25fcc 
(green shaded). 

Our simulations show that the cortical contributions to sur- 
round suppression and contrast dependent receptive field size 
account for a large fraction, though not all of the magnitude 
of the phenomena observed experimentally. This leaves room 
for contributions from the LGN input. It seems reasonable 
to assume that LGN cells in macaque will display both extr- 
aclassical phenomena, although, somewhat surprisingly, this 
has to our knowledge not been verified yet. Surround sup- 
pression of LGN cells at cortical optimal spatial frequencies 
has been observed in marmoset 16 and cat 9 . Contrast depen- 
dent receptive field growth of LGN cells has been observed in 
marmoset and an average growth ratio of 1.3 was reported 16 . 
Transfer of receptive field growth of LGN cells to VI corti- 
cal cells seems unavoidable because it simply introduces an 



overall scaling factor on the entire visual input in VI. In- 
heritance of LGN surround suppression by VI cells is not so 
obvious. Consistent with the Hubel and Wiesel view (the way 
our model is constructed), LGN input arrives in a VI cell es- 
sentially as output of small clusters of about 10-20 LGN cells. 
It is thus not immediately clear if and how much LGN sur- 
round suppression can be transferred to VI cortical cells. We 
addressed this issue by repeating our simulation, originally 
performed at kc, at the smaller spatial frequency of 0.25fcc- 
Results of this simulation are shown in in the remainder 
of Figure 8. We see (Fig. 8B) that practically all LGN sur- 
round suppression is transferred to cortical cells. Given that 
an average surround suppression of SI\ ~ 0.4 is observed in 
macaque 8 , our results show that cortical short-range con- 
nections together with surround suppression present in LGN 
cells can easily explain the degree of surround suppression 
seen experimentally. We also see that the presence of LGN 
surround suppression also has consequences for the prevalence 
of the different mechanisms by which cortical suppression is 
achieved. When LGN surround suppression is present, the 
prevalence of the suppression mechanisms is substantially al- 
tered in favor of mechanisms B and C (which require a reduc- 
tion of excitation) at the expense of mechanism A (increased 
direct inhibition). Results for contrast dependent receptive 
field size are practically unaltered when LGN suppression is 
included (not shown). 

The cortical magnification factor 

As mentioned in the Introduction, arguments in favor of in- 
volvement of long-range connections and/or extrastriate feed- 
back in extraclassical phenomena, are indirect and all rely on 
the cortical magnification factor as a key ingredient. Recep- 
tive field size and scatter are systematically ignored. It is ar- 
gued that surround sizes would be too big to result from local 
short-range connections. We have already shown that, on the 
contrary, through polysynaptic interactions in the network it 
is possible to create the surround sizes observed experimen- 
tally with only local short-range connections. 

One naturally wonders to what extent our findings depend 
on the actual value of the cortical magnification factor. Intu- 
itively, a smaller cortical magnification factor is not beneficial 
for the role of short-range connections in the creation of extr- 
aclassical receptive field phenomena, since these connections 
cover less visual space. However, the minimum amount of vi- 
sual space covered is set by the receptive field size and scatter. 

To check whether this minimum visual range of cortical 
short-range connections is in itself sufficient to generate sur- 
round suppression and contrast dependent receptive field size 
we repeated our simulations with an infinite cortical magnifi- 
cation factor, v~ l = oo (geometric parameter Q = oo, all else 
unchanged, see Methods,). The results are shown in Figure 9 
for the P0 case (M0 yields similar results). Clearly, the finite 
receptive field scatter by itself, 60% (M0) and 30% (P0) of the 
average receptive field size (Methods), is sufficient to generate 
both extraclassical phenomena to practically the same degree 
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as it does in the presence of a realistic cortical magnification 
factor. Other properties of our model discussed in this paper 
are also not qualitatively different for infinite cortical magnifi- 
cation factor, with one exception. As is apparent from Figure 
9A, it is now more difficult to connect the LGN axons in such 
a manner that the organization of orientation preference and 
ocular dominance displays the same level of order as seen for 
a realistic cortical magnification factor (Fig. 1). (This does 
of course not imply that order could not be improved with 
more specific connections than the ones used). 

Arguments ruling out local short-range connections and 
LGN input as the origins of extraclassical phenomena based 
on the cortical magnification factor are inherently weak, since 
it is a macroscopic measure and thus inferences based on it 
regarding which cells could influence which cell and on what 
time scale, cannot be very precise. Our simulations in this 
section further challenge such arguments, by showing that re- 
ceptive field size and scatter by themselves, regardless of the 
cortical magnification factor, can be a determinative factor 
for extraclassical receptive field phenomena. 




Figure 9: Results for infinite cortical magnification factor (fi = 
oo ) for the PO configuration, MO configuration yields similar re- 
sults, (left) Simulated optical image of ocular dominance and ori- 
entation preference, in the spirit of optical imaging experiments 28 
(compare Figure 1). (right top) Suppression index 5*11 at high (un- 
filled) and low (green shaded) contrast (compare Figure 2). (right 
bottom) Receptive field and surround growth ratios, r-/r+ (blue 
shaded) and R-/R+ (unfilled). 



Discussion 

There is considerable debate over the origins of extraclassical 
receptive field phenomena such as surround suppression and 
contrast dependent receptive field size. A primary reason is 
that there are no experimental data yet available that directly 
point to mechanisms of these phenomena. This requires that 
any data needs to be interpreted through some sort of model 
or theoretical framework. At the same time, such a model 
or theory needs to be sufficiently sophisticated, for example, 



it needs to adequately address classical response properties 
since they set the context for the extraclassical phenomena. 
We believe our model is the first anatomically and physiolog- 
ically realistic model to simultaneously address classical and 
extraclassical spatial summation. 

We know of only one attempt in the published literature 35 
to develop a spiking neural network model to address the extr- 
aclassical response phenomena discussed in this paper. That 
model does not address classical responses, and neural mech- 
anisms for extraclassical phenomena, by construction, arise 
from long-range connections. Further, in that model contrast 
dependent receptive field size is achieved via contrast depen- 
dent surround suppression, that is, the surround suppression 
systematically decreases as a function of decreasing contrast, 
which contradicts experimental data 4 ' 8 . 

We have shown that considerable surround suppression 
and contrast dependent receptive field growth can be sponta- 
neously generated solely by short-range cortical connections in 
VI without any contributions from other sources. We demon- 
strated that surround suppression of LGN cells, if present, is 
easily transferred to VI. Our simulations thus provide rigor- 
ous computational support for the intriguing hypothesis that 
the LGN input and cortical short-range connections in VI 
are primarily responsible for the phenomena, with little or 
no contributions from long-range or extrastriate connections. 
We showed that with only 25% suppression in the LGN cells, 
the cortical surround suppression in the model exceeds the 
suppression observed experimentally. More radically inter- 
preted, our results thus suggest that long-range lateral con- 
nections and/or extrastriate feedback contribute negatively to 
surround suppression, that is, rather than being suppressive, 
our results suggest that their contributions are in fact facili- 
tatory. 

In our model all three neural mechanisms for surround sup- 
pression are active, in agreement with experimental observa- 
tions in cat 7 . When LGN suppression is included, we observe 
strong contributions to cortical suppression from a reduction 
of recurrent cortical excitation, rather than from an increase 
in direct cortical inhibition. We find, on average, a growth 
of spatial summation extent of excitation and inhibition at 
low contrast, as predicted by DOG 4 and ROG models. But 
this growth bears no simple relationship with the receptive 
field growth seen in spike responses, which usually involve 
other/additional changes in the relative gain of the excita- 
tory and inhibitory inputs. Notably, significant receptive field 
growth is usually much larger than the growth of the spatial 
summation extent of excitation and inhibition (Fig. 7). 

As does the biological primary visual cortex, our model 
produces these properties in distinctly different geometric set- 
tings (with identical strength parameters), namely for the 
magno and parvo input layers, at para-foveal eccentricities 
and around 10° eccentricity. Dimensional observations (pa- 
rameter Vt) imply our model results also translate to 30° ec- 
centricity, modulo a geometric scaling factor. The ubiqui- 
tous nature of these phenomena and their mechanisms in our 
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model suggest that they are basic response properties of VI. 
Given realistic classical response properties, as discussed in 
this paper, they seem to require little more than receptive 
field scatter and isotropic short-range connectivity, with per- 
haps some weak constraints on macroscopic organization of 
ocular dominance and orientation preference. They do not 
seem to require more elaborate architecture or physiological 
properties, such as specific cortical connectivity, long-range 
connections within VI, extrastriate feedback, synaptic depres- 
sion/facilitation etc.. 

Recent data for cat suggest that partial inheritance of sur- 
round suppression in VI from LGN cells does indeed occur 9 . 
Other recent experiments show that strong surround suppres- 
sion is observed for drifting gratings having spatial and tem- 
poral frequencies outside the range at which most cortical 
cells typically respond, indicating these signals arise within 
the input layer of VI or the LGN itself 36 . Given the general 
nature of our results, we may conclude that the presence of 
the phenomena in LGN cells could in principle be of the same 
origin as we have suggested here for VI cells. 

There are various further aspects of surround suppression 
that we have not explicitly addressed in this paper. Among 
them are orientation tuning of the surround and dynamics 
(timing) of the suppression. Preliminary simulations indicate 
that orientation tuning of the surround is well-captured by 
our model. For what concerns timing of the suppression seen 
in our model, it is clear that through polysynaptic interac- 
tions in the network, delays of the onset of surround sup- 
pression could range anywhere from 0-20 ms, conservatively 
estimated. Indeed, preliminary simulations show that timing 
of the surround suppression in our model is consistent with re- 
cent experimental findings. Our model has a rich dynamics 
and is well-suited to also yield relevant results regarding the 
dynamics of surround suppression. This is one of our interests 
for future research. 
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Supplementary Information 

Model cortex 

Our model cortex consists of N = 256 x 256 conductance 
based integrate-and-fire point neurons (one compartment), of 
which 75% are excitatory and 25% are inhibitory and which 
are randomly distributed on a square grid. 

Rcscaling of Eq. (1) in the paper works as follows. Nu- 
merical values for the biophysical parameters are, for the 
capacitance C = 10 -6 F cm~ 2 , the leakage conductance 
gL = 50 x 10~ 6 fl^ 1 cm -2 , the leakage reversal potential 
Vl = —70 mV, the excitatory reversal potential ve = mV, 
the inhibitory reversal potential vj = —80 mV, the spiking 
threshold vt = —55 mV, and the reset potential vr = —70 
mV. The equation is then divided by C *Vt- Physiological 
conductances (SI -1 cm~ 2 ) and currents (Ampere cm~ 2 ) are 
obtained by multiplication with 10~ 6 and 15 x 10~ 9 respec- 
tively. 

The functions C ,' (r) in Eq. (3) of the paper describe the 
cortical spatial couplings (cortical connections), they are con- 
structed as follows. We assume the availability of postsynap- 
tic sites Nd on a cell (dendrites) to decay exponentially as a 
function of distance with length scale D^,' 1 ^, that is, Nd <~ 
cxp[— (r/D^,' k ^) 2 ], and make a similar assumption for the 

presynaptic sites A^ a (axons), N a <~ exp[— (r/A ,' ) 2 ]. Then 
the spatial coupling strength (assuming individual synapses 
have equal strength) between two cells decays exponentially 
with length scale (a^) 2 = (£#£) 2 + {A*,*) 2 and can be 
written as 

<::m = v;X : ;> x p[-^/v:") 2 ]' ( 17 ) 

with the normalization constants 

Ki = l E cx pHINI/v;") 2 ]} • ( 18 ) 

k' k 

In this way, the parameters c ,' are interaction strengths that 
define the density and length scale invariant contribution of 
population Vkip) to the conductance of a cell in population 
Vy (m')- The change in membrane potential of cell i € Vk 1 (aO 
due to a single spike of cell j e Vkip) is proportional to 
V^V^^'^^^^'exph^/aJ;^) 2 ], where n k ^ is the 
cell density of population P k (n) and r it j = \\xi — Xj\\. 

The temporal kernels G MJ (r) in Eq. (3) of the paper de- 
scribe the synaptic dynamics of cortical synapses, are normal- 
ized to unity, G^.j(T)dT — 1, and arc of the form 

GW(r) = 
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if T < 

&am (re- T / a ">') 5 if < t < A^a^ . 

(A M a M . ie - A -) 5 e -<- A .<w)/i>» if T > A^a^j 

(19) 

The time constants are based on experimental observations 38 . 
The excitatory kernels Gg,., have a fast (AMPA) component 
defined by the peak times a,E,ii which are drawn from a uni- 
form distribution between 1 ms and 4 ms, and a slow (NMDA) 
component which is defined by the decay time 6g = 15 ms. 
Transition between the two regimes is set by A^ = 4/3. Sim- 
ilarly, the inhibitory kernels G/,i have a fast (GABA) compo- 
nent set by a^i, chosen from a uniform distribution between 
3 ms and 6 ms, and a slow component 39 denned by 6/ = 10 
ms, while A/ = 3/2. The constants k^i are normalization 
constants. These kernels imply a spike memory on the order 
of 70 ms for excitation and 50 ms for inhibition. 

The cortical spatial coupling length scales are taken to be 
in agreement with anatomical data 40-43 , axon and dendrite 
parameters are A ,' E = 200 fim (axons of excitatory neu- 
rons), A j = 100 fim (axons of inhibitory neurons), and 

D k ,' k — 50 fj/m (dendrites of all neurons) respectively. This 
implies that excitatory connections, both on excitatory and 
inhibitory cells, have characteristic length scale a ,' E ~ 200 
/i 77i while inhibitory connections, both on excitatory and in- 
hibitory cells, have a characteristic length scale a ~ 100 
fim. So, in agreement with the anatomy our model has in- 
hibitory connectivity which is of a much shorter spatial range 
(a factor of two) than the excitatory connectivity. 

Because of lack of experimental data, strength parameters 
in the model are free parameters. They occur in the LGN 
input, noise levels, and the cortical interaction. In particular, 
the cortical coupling strength matrix c ,' consists entirely of 
such parameters. It is currently not experimentally feasible 
to directly determine the components of the cortical coupling 
strength matrix c ,' or the other strength parameters in the 
model. A natural question is "which subsets of strength pa- 
rameters from the set of all possible strength parameters yield 
physiologically realistic results?" . A complete answer to the 
question is also not yet feasible, because of the size of the 
simulation and the number of parameters involved. 

We provide in this paper a set of strength parameters that 
yields physiologically realistic results. The set of strength pa- 
rameters we provide is obtained by adhering to a few general 
principles regarding the LGN input and cortico-cortical in- 
teractions. These general principles are: (a) No distinction 
between LGN input in the excitatory and inhibitory cell pop- 
ulations, (b) Cells with LGN input, both excitatory and in- 
hibitory cells, receive their excitation in about equal amounts 
from LGN input, cells with LGN input, and cells without 
LGN input, (c) Cells without LGN input, both excitatory 
and inhibitory cells, receive most of their excitation from cells 
with LGN input, (d) Cells with (without) LGN input, both 



excitatory and inhibitory cells, receive most of their inhibition 
from cells with (without) LGN input. 

Proceeding in this way, we considered a small number of 
classical response properties for setting the operating point 
(set of strength parameters). Specifically, these classical re- 
sponse properties are: (i) Absence of any global phase-locked 
oscillations and synchrony, both under visual stimulation and 
without visual stimulation. This requirement limits the over- 
all maximum size of the strength parameters, (ii) Distribu- 
tion of activity (firing rates) over the cell population, with and 
without visual stimulation. This requirement constrains the 
overall balance between excitatory and inhibitory strength pa- 
rameters, (iii) Distribution (over cell population) of response 
modulations for drifting grating stimuli. This essentially is 
a requirement on the model's composition in terms of sim- 
ple and complex cells, (iv) Orientation tuning of both simple 
and complex cells. This requirement results in large conduc- 
tances, particularly, large inhibitory conductances. Note also 
that this requirement is a major constraint, since orientation 
tuning for drifting gratings is not easily achieved, e.g. see 18 . 
This is so because of the fact that the average LGN input is 
practically untuned for orientation and inhibitory connections 
are of much shorter range spatially than excitatory connec- 
tions. 

Numerical values for the strength parameters were obtained 
by performing many trial-and-error simulations, in a properly 
scaled-down version of the model, in search of a suitable oper- 
ating point. The coupling matrix resulting from this approach 
and used in the numerical simulations presented in the paper 
is 
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This coupling strength matrix (Eq. 1201) and the other 
strength parameters (in noise levels and LGN input, to be 
discussed below), were chosen because they set the operat- 
ing point of the model such that it reproduces the classical 
response properties (i)-(iv) in agreement with experimental 
data. Further support for the physiological relevance of this 
operating point is provided by the fact that, in this setting, 
model properties accurately extrapolate to a variety of other 
known properties of the biological visual cortex beyond prop- 
erties (i)-(iv). Among them are (v) spatial and temporal fre- 
quency tuning (vi) a cortex operating at high conductance 
levels, (vii) distributions of response modulations for contrast 
reversal gratings, (viii) dynamics of responses to stochastic 
stimuli (reverse correlation), and, of course, (ix) the extra- 
classical response properties which are the main topic of the 
paper. Furthermore, results obtained at this operating point 
are robust, that is, they remain quantitatively similar and 
thus physiologically realistic, when a random change < 10% 
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is introduced in the set of strength parameters. 

The external stochastic terms rj^^it) in Eq. (1) of the paper 
are given by 

/oo 
G^(t-r)S^(T)dr . (21) 
-oo 

Where the kernels have the same form as (|19f) (with 
randomly selected a^'s) and are Poisson spike trains 
(mean firing rates 100 spikes/s (/i = E) and 125 spikes/s 
(p = /)) belonging to neuron i (different ones for each cell). 
The noise strengths 77° i are drawn from a uniform distribution 
between 1 and 5 if i £ Vo(E), are equal to 2 if i £ Vi(E), 
are drawn from a uniform distribution between and 30 if 
i £ Vq(I), and are drawn from a uniform distribution between 
16 and 46 if i £ Vi(I). 

Geometric parameters and LGN input 

The sets Nq G - n are constructed as follows. Our 4x4 mm 2 
modelcortex is partitioned into 8 parallel bands (0.5x4 mm 2 ), 
which alternate representation between the two eyes. Subse- 
quently, initial retinotopic maps for each eye were defined as 
the identity map plus scatter as follows. The 30% (members of 
Vi(/j,)) of the neurons that receive LGN input, were assigned 
initial RF centers (yk) depending on their positions (xk) in 
the model cortex via yk — vxk + jpk, where v is the inverse 
cortical magnification factor, 7 is the initial RF scatter and 
pk are scatter vectors, components of which are drawn from a 
uniform distribution on [—1, 1]. The parameter v is 0.2° /mm 
for 0° eccentricity (M0,P0) and 0.7° /mm for 10° eccentricity 
(M10,P10), which correspond to cortical magnification fac- 
tors at the lower end of their experimentally observed range 
in macaque 44 ~ 46 . The scatter parameter 7 is 0.3°, 0.1°, 0.35° 
and 0.35° for M0, P0, M10 and P10 respectively. These scat- 
ter values are in the experimentally observed range 44 ' 47 and 
moreover assure a more or less uniform initial distribution of 
cortical receptive fields. 

Next, the neurons that were assigned initial receptive field 
centers are connected to LGN axons of the corresponding 
eye to set up their orientation preference (pinwhccls). Be- 
sides an initial RF center, each cell in Vi(p) (for each eye) 
is assigned a template for the organization of the ON and 
OFF sub fields of its initial RF, which is randomly chosen 
from the 4 basic symmetry configurations seen experimen- 
tally 47 . LGN-cortical connections are initially made so as 
to best approximate the assigned initial receptive field center 
and template for each cell. Then, the LGN-cortical connec- 
tions are rearranged by disconnecting and reconnecting cells, 
so as to achieve LGN axon sizes that agree with the anatomi- 
cal findings for macaque 29 ' 30 . For each disconnection, a new 
connection is selected from all possible LGN axon candidates 
(in cortex), such that the LGN cell belonging to this newly 
selected axon has a receptive field closest to the receptive field 
of the LGN cell that was disconnected. If no candidates were 
available no reconnection was made. The constraint put on 



the sizes of the LGN axonal arbors in agreement with their 
experimentally observed maximum sizes is important, since 
we find that this puts considerable restrictions on the possi- 
ble connections, leading to a realistic local scatter in preferred 
angle (as observed experimentally 48 ). 

In our description of the LGN input we use the term "ini- 
tial" when referring to receptive fields and retinotopic maps, 
since these were only assigned to build the connections be- 
tween LGN cells and a fraction (30%) of the cortical cells. 
Actual receptive fields and retinotopic maps in our simula- 
tions depend on the final connections. In addition, they are 
largely shaped by the cortical interactions, since there are 
many cells (70%) that do not receive LGN input. The actual 
retinotopic maps are approximately identity maps, but with 
scatter somewhat smaller than in the initial maps. We see a 
rich variety in shapes, sizes and organization of ON and OFF 
sub fields, much like what is seen experimentally 47 . The re- 
ceptive field sizes and spatial and temporal frequency tuning 
properties of our model cells also agree well with experimental 
data. 

Data collection 

Preferred angles, spatial and temporal frequencies are ob- 
tained by using a drifting grating stimulus for one eye (other 
eye 7 = 0), with high contrast (e = 1), in a large aperture 
(ta ~ 7-10 times the average receptive field size). Responses 
to contrast-reversal stimuli, i.e. cos( u>t — k ■ y + <fi) of the 
drifting grating stimulus is replaced by cos( cut) cos(fc • y + cf>), 
are obtained in a similar way. The receptive field centers are 
mapped using a drifting grating stimulus for one eye (other 
eye 1 = 0), with high contrast (e = 1), in a small aperture 
(r\4 ~ l/4th of the average receptive field size). The apertures 
are centered on a rectangular grid (5x5 or 6x6, grid spacing 
about half of the average RF size) which itself is centered 
on the visual field covered by our model cortex. The angle, 
spatial and temporal frequencies of the grating are kept fixed 
during this experiment. The temporal frequencies are set to 
the averaged preferred values for the case under consideration. 
The preferred parameters and the receptive field centers are 
from averaged spike responses, using the mean response for 
complex cells and the first harmonic for simple cells. Each 
stimulus presentation (all parameters fixed except time) lasts 
for 4 s and is preceded by a 1 s blank stimulus (e = 0). This 
procedure is repeated twice, starting with different initial con- 
ditions and different external noise realizations. Average re- 
sponses are obtained by averaging over cycles and trials. 

The primary data, i.e. responses and conductances as a 
function of aperture size for single eye stimulation, are ob- 
tained with the temporal and spatial frequencies of the grat- 
ing set equal to the averaged preferred values for each case 
(M0, M10, P0, P10). Data samples consist of cells that have 
their preferred angle equal to the grating angle, their preferred 
temporal frequency within 2 Hz of the grating frequency, a 
preferred spatial frequency k p that satisfies < k v < 2k 
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where k is the grating spatial frequency, a receptive field cen- 
ter that is less than 4jth of the average RF size away from the 
aperture center, a maximum response at low contrast that is 
greater than /& + 5 where ff, is the mean blank response (in 
spikes/s), and finally, a central cortical location confined to 
the dashed white rectangle in Figure 1 of the paper. In this 
way we collect approximately 200 cells in each sample. 

The experiments were performed at "high" contrast, e = 1, 
and "low" contrast, e = 0.3. In our model contrast depen- 
dence of the LGN cells is approximately linear. Most cortical 
model cells, consequently, also have a fairly linear contrast 
response function. Therefore, the numerical contrast values 
used in the model should not directly be compared with nu- 
merical contrast values used in experiments with animals, but 
rather they should be compared only with the linear part of 
the contrast response function, with e = 1 being the "high" 
part and e = 0.3 being the "low" part. Note that this also 
means complications resulting from response saturation at 
high contrast are not present in our model. Indeed, our find- 
ings presented in this paper are insensitive to using instead of 
e high = 1 a lower value of e n jgh = 0.8, for which high contrast 
responses are approximately 80% of their saturation values. 
Further, when we model the contrast response functions of M- 
LGN cells more realistically, using experimental data from 16 , 
and confine our experiments to the linear part of the contrast 
response functions, our findings remain unchanged. 

Simplified model equations 

We rewrite Eq. (1) of the paper as 19 



dt 



-g T ,k(t,r A ) v k + I D ,k(t, r A ) 



(22) 



where 

9T,k(t,r A ) = 9l + 9E,k{t, [S]e,ta) + 9i,k{t, [S]i,r A ) (23) 

lD,k(r A ,t) = 9E,k{t, [S] E ,r A ) V E - gi,k(t, [S]i,r A ) \Vi\ . 

(24) 

Recall that our model is in fact a system of many strongly 
coupled equations. Together they determine the total con- 
ductance gT,k(t,r A ) and difference current lD,k(t,r A ) which 
govern the behavior of a model cell. Our model cortex oper- 
ates at large conductance levels as part of the requirement to 
achieve a broad range of known properties of the biological 
visual cortex. Indeed, several experiments have demonstrated 
large cortical conductances under in vivo visual stimulation, 
see e.g. 49 for a review. A large total conductance implies 
that the membrane potential is well-approximated simply by 
Ohm's law v k (t,r A ) w I Dyk (t,r A )/g T . k (t,r A ). We also find 
in our model that the cycle-to-cycle and trial-to-trial fluctu- 
ations in gx,k(t,r A ) are considerably smaller than the mean. 
Therefore, we find that to a good approximation 19 



(v k (t,r A )) w V k (r A ,t) = -- 

{9T,k{t> r A)) 



Further, we find that for our model cells the spike rates are 
well-predicted by a rectification model which assumes the fir- 
ing rate to be proportional to the rectified total current at the 
spike threshold (v k = 1), 

(S k (t,r A )) fk(t,r A ) = 



S k [{I D ,k{t,r A )) - (0T,H(t,r A )) ~ A fe ] + , (26) 

where [x] + = x if x > and [x]+ = if x < 0, and where, to 
a good approximation, the gain S k and threshold A k do not 
depend on the aperture radius r A for most cells. With use of 
expression <|25|) for the membrane potential, our rectification 
model (Eq. I26fl is equivalent to the standard rectification 
model 50 (S k (t,r A )) « S k [(v k (t,r A )) - v T . k } + , where the 
gain parameter 5 k and threshold vx,k however depend on the 
total conductance, that is, on aperture size and time. The 
standard rectification model with constant (cell specific) gain 
and threshold parameters has been shown to work well in cat 
51 for large drifting grating stimuli of varying angles. Our rec- 
tification model gives a significantly better fit for our stimuli 
than the standard rectification model (with fixed gain and 
threshold). Some fits are shown in Figure 10. It provides 
a relation between firing rate and excitatory and inhibitory 
conductances in which these conductances act simply via ad- 
dition and subtraction, respectively. Since the parameters are 
constants it allows us to base our analysis of the firing rate 
(as with the membrane potential) directly on the behavior of 
the excitatory and inhibitory conductances as a function of 
aperture size and time. 

Simulated optical imaging 

To simulate optical imaging experiments 28 for ocular domi- 
nance and orientation preference, high contrast drifting grat- 
ings in a large aperture are presented to the left and right eye 
separately, temporal and spatial frequencies of the grating are 
set equal to averaged preferred values. Stimulus presentation 
is as discussed for finding the preferred angles etc.. We as- 
sume the optical signal at pixel £ arising from a single eye (L 
or R) stimulation with a grating with angle 9 k to be 
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Q = L or R , (27) 



where Tj^Q(0 k ) is the mean (spike) response of cell j and N n (t) 
is the neighborhood of the nearest n cells to pixel I. Using this 
assumption we can produce an arbitrary "optical" imaging 

result for our model by further processing of t^q (9 k ; n) in the 
same spirit as in optical imaging experiments. We compute 
ocular dominance Oi{n) for a given pixel of our model cortex 
as 



Oi(n) 



(25) 
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Figure 10: Some examples of the approximations given by Eq. 
(6) and (7) of the paper for a simple cell (left) and a complex 
cell (right) for apertures around the aperture of the maximum 
response. Figure layout is as Figure 6 of the paper. (A & B) 
Comparison of the instantaneous firing rate (<S(i, ta)) to the rec- 
tification model approximation ,f(t,rA) (Eq. 1261 . Note that the 
approximation is reasonable for all aperture sizes, but better for 
those that impinge on the suppressive surround. Fitting parame- 
ters are S = 1.2, A = -100 s _1 (left) and 5 = 1.8, A = -25 s _1 
(right). (C & D) Comparison of model responses {v(t, Ta)) with 
the Ohmic (slaving) approximation V(t,rA) (Eq. 12 "it . indicating 
the model's operation at high conductance states. 



Orientation preference Qp.i{n) (preferred angle of pixel £) is 
computed from the "averaged difference vector" similarly as 
in actual experiments 28 , 



? (n) = -arg^i ^ i 

k (Q=L,R 



(29) 



Results in Figures 1 and 9 of the paper are obtained with 
ri = 75. 

Extraclassical spatial summation 

A summary of our data for the four configurations MO, M10, 
P0, and P10 is provided in Figure 11. The figure is organized 
in columns and rows, each column corresponding to a partic- 
ular configuration, and each row corresponding to a particu- 
lar response measure. All results shown are based on spike 
train responses. The second row shows the distributions of 
receptive field size for high and low contrast (r + ._), the third 
row shows the distribution of the surround size at both con- 
trasts (i? +i _). The distributions of receptive field size and 
surround size show good agreement with experimental data 
8 . Notice the growth of both receptive field size and surround 
size for low contrast. Also, note the considerable diversity in 
receptive field and surround sizes in the model. This is not a 
result of any variability in the LGN receptive field sizes, all 



LGN cells for a particular configuration have identical recep- 
tive field size. Instead, the diversity seen in the receptive field 
and surround sizes of the cortical cells is a result of the corti- 
cal interactions and realistic constraints on LGN axon sizes. 
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Figure 11: Summary of extraclassical spatial summation for the 
four model configuration. In rows 1-3 unshaded histograms are 
for gratings at high contrast (+), shaded at low contrast ( — ). 
Throughout the figure, means are indicated by arrows, solid ar- 
rows refer to shaded histograms, dashed arrows to unfilled his- 
tograms. (ROW 1) Distributions of the suppression index SIi. 
All suppression is exclusively due to short-range cortical connec- 
tivity. (ROW 2) Receptive field size distributions (r+__). Dis- 
tributions are shifted to the right for low contrast indicating an 
increase in receptive field size. Receptive field sizes in the model 
are in agreement with experimental data 8 . (ROW 3) Surround 
size distributions (_/?+,_) for high and low contrast. Surround sizes 
are in agreement with experimental data 8 . (ROW 4) Histograms 
of the ratios of the receptive field and surround sizes at low and 
high contrast, r_/r+ (blue shaded) and R-/R+ (unfilled). Growth 
of surround sizes is about the same as for receptive field sizes and 
each case shows a similar, significant growth of about a factor of 
two (Wilcoxon test on ratio larger than unity: p < 0.001 for both 
receptive field and surround growth). 



Mechanisms of contrast dependent receptive 
field size 

Consider a situation where both gE and gj have their maxima 
at the same aperture size tb = fj = r* and are monotonically 
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increasing for ta < t* and monotonically decreasing for ta > 
r*, as depicted in Figure 12A. We can distinguish three classes 
with respect to the relative location of the maxima in spike 
responses r$ and the conductances r*, namely {X: 7^5 < r*}, 
{Y: rs = r*} and {Z: r$ > r*}. It follows from that 
if we define the parameter Gq(v) = (\vj\ + v)/(ve — v) then 
we can characterize the difference between classes X and Z by 
the way that G crosses Go(l) around rs as depicted in Figure 
12 A. For class Y the parameter G is not of any particular use 
as it can assume arbitrary behavior around rg. It follows from 
(I12fl that similar observations hold for the maximum in the 
membrane potential r v and we need simply to replace Gq(1) 
with Go(v(r v )). Obviously, a growth of receptive field can 
occur without any change in the spatial summation extent 
(r*) of the conductances. Suppose we wish to remain within 
the same class X or Z, then receptive field growth, can be 
induced, for instance, by an overall increase (X) or an overall 
decrease (Z) in relative gain G{ta) as shown in Figure 12 A 
(dashed line). Receptive field growth also can be caused by 
more drastic changes in G so that the transitions X — * Y, 
X — > Z or Y — * Z occur for a high — > low contrast change. 
The situation is somewhat more involved when we allow for 
non-suppressed responses and conductances, and for different 
positions of the maxima of qe and gi, however, the essence 
of our conclusions remains the same. 

From a more precise analysis based on the relative gain 
parameter G, we find that for more than 50% of the cells with 
significant receptive field growth, a transition takes place from 
a high contrast RF size less or equal to the spatial summation 
extent of excitation and inhibition, to a low contrast receptive 
field size which exceeds both. This analysis is summarized in 
the remainder of Figure 12. 

Cells were classified (Fig. 12B) according to the relative 
positions of their maxima in spike response (rs) and exci- 
tatory (r E ) and inhibitory (77) conductances, using F0+F1 
components. Membrane potential responses yield similar re- 
sults. Comparing this classification at high and low contrast 
we observe a striking difference for cells with significant re- 
ceptive field growths, i.e. with growth ratios >1.5 (Fig. 12B, 
bottom), indicative of X — > Y, X — > Z and Y — > Z transitions 
(as discussed in the simplified example above). 

In this realistic situation there are of course many more 
transitions (i.e. 13 2 ), however, that we indeed observe a 
prevalence for qualitatively these transitions can be demon- 
strated in two ways using slightly modified definitions of the 
X,Y,Z classes. First (Figure 12C, left), if we redefine the 
X,Y,Z classes with respect to rs and 775 while ignoring 77, 
i.e. {X: rs < 77;}, {Y: rs = rs} and {Z: rs > rs}, then the 
transition distribution for cells with significant receptive field 
growth shows that in about 60% of these cells a X — > Z or 
Y — ► Z transition occurs. Taken together with the fact that 
roughly 10% of the cells with significant receptive field growth 
(Figure 12B, bottom) have 77 < rs < rE at high contrast and 
te < rs < 77 at low contrast, we can conclude that for more 
than 50% of the cells with significant receptive field growth, 



a transition takes place from a high contrast RF size less or 
equal to the spatial summation extent of excitation and in- 
hibition, to a low contrast receptive field size which exceeds 
both (by at least one aperture). Note that these transitions 
occur in addition to any growth in or r/. Secondly (Figure 
12C, right), the same conclusion is reached when we redefine 
the X,Y,Z classes with respect to rs and 77 while ignoring 775 
({X: r s < 77}, {Y: r s = 77} and {Z: r s > 77}), Now aX-t 
Z or Y — > Z transition occurs in about 70% of the cells with 
significant receptive field growth, while about 20% of the cells 
with significant receptive field growth (Figure 12B, bottom) 
have 775 < rs < rj at high contrast and 77 < rs < rE at low 
contrast. 

Figure 12C also demonstrates the presence of a rich diver- 
sity in relative gain changes in our model, since all transitions 
(for all cells, unfilled histograms) occur with some reasonable 
probability. Finally, Figure 12C establishes that there is a 
relationship between the difference in the prevalence of the 
surround suppression mechanisms at high and low contrast 
and receptive field growth. To see this, first note that for 
the redefined Y and Z classes with respect to rs and 77, the 
surround suppression, if any, must be caused by mechanisms 
B or C. Thus, since the total probability (Figure 12C, right) 
of transitions • — * Y and • — ► Z (not including Y — * Y, Y 
— * Z and Z — > Z) is clearly larger than for transitions • — * X 
(not including X — ► X), this means that surround suppression 
mechanisms B and C must be more prevalent at low contrast 
than at high contrast, which is confirmed by our data (not 
shown) . 
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Figure 12: (A) Schematic illustration of mechanisms for receptive 
field growth under equal and constant spatial summation extent 
of the conductances (t\e = Ti = r*). (B) Distributions of the 
relative positions of the maxima (receptive field sizes) of spike re- 
sponses rs and conductances te and rj, for the MO configuration 
(other cases give similar results). A division is made with respect 
to the maxima in the conductances, this corresponds to the left 
(te = ri), central (te > ri), and right (te < rj) part of the 
figure. Each panel is further subdivided with respect to the max- 
imum in the spike response rs- Upper histograms are for all cells 
in the sample, lower histograms are for cells that have receptive 
field growth r_/r+ > 1.5. Unfilled histograms are for high con- 
trast, shaded histograms are for low contrast. (C) Prevalence of 
transitions between positions of maxima in spike responses and ex- 
citatory conductances (left) and in spike responses and inhibitory 
conductances (right) for a high — > low contrast change. See text 
for definitions of X, Y, Z classes. Data are evaluated for all cells 
(unfilled histograms) and See text for for cells with a receptive field 
growth r_/r+ > 1.5 (shaded histograms). 
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